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As two important mechanical properties, strength and ductility generally tend to be muturally exclusive in 
conventional engineering materials. The breakthrough of such a trade-off has been potentiated by the 
recently developed CuZr-based bulk metallic glass (BMG) composites ductilized by a shape memory 
CuZr(B2) phase. Here the microstructural dependences of tensile properties for the CuZr-based BMG 
composites were elucidated qualitatively and modeled quantitatively, and the underlying mechanisms were 
unraveled. Through the microstructural percolation induced by matching the length scales of particle size 
and interparticle spacing, a notable breakthrough was achieved in the composites that the general conflicts 
between strength and ductility can be defeated. This study is expected to greatly aid in the microstructural 
design and tailoring for improved properties of BMG composites. It also has implications for the 
development of strong and ductile materials in the future. 

As two important mechanical properties, strength and ductility usually tend to be mutually exclusive in 
many classes of materials with the general rule that they are inversely related 1 " 3 . The changes in a 
material's microstructural architecture often affect the strength and ductility in opposite ways, i.e. 
strengthening but sacrificing the ductility or alternatively ductilizing at a cost of strength. Thus a general 
trade-off between the two properties seems inevitable to make a compromise in developing strong and ductile 
materials. In particular, so far the structural applications of explored high-strength materials are usually limited 
by their relatively low ductility and toughness in engineering practice 3 . This is also the case for the emerging bulk 
metallic glasses (BMGs) which usually demonstrate extraordinarily high strength but disappointingly low duct- 
ility 4 ' 5 . The wide -spread structural applications of BMGs have been strictly limited by their macroscopic room- 
temperature brittleness and strain softening nature. It is revealed that these shortcomings can be mitigated by 
combining the glassy matrix with in situ formed ductile crystalline phases, i.e. forming BMG composites 6 ' 7 . In this 
respect, the CuZr-based BMG composites ductilized by a shape memory CuZr(B2) phase have attracted great 
attention in recent years owing to their outstanding mechanical properties 8 " 12 . The unique transformation- 
mediated deformation mechanisms endow these composites with pronounced tensile ductility and work-hard- 
ening capacity without significantly scarifying the high strength of the glassy matrix 8 ' 9 . In this scenario, the 
successful combination of shape memory alloys with BMGs may provide new opportunities to defeat the general 
conflicts between strength and ductility and open new avenues for the development of high-performance struc- 
tural materials. 

It is known that the mechanical properties of BMG composites depend strongly on their microstructures 10-13 . 
Thus the elucidation of the structure-property correlations is invariably essential for the microstructural design 
and tailoring to achieve intended mechanical properties. However, hitherto such an issue has rarely been dis- 
cussed and still remains far from understood in terms of the tensile properties for the CuZr-based BMG 
composites 11 . As a consequence, it still remains unclear and challenging about how to achieve the breakthrough 
of the trade-off relation between strength and ductility in the composites through controlling their microstruc- 
tures. On the other hand, the percolation theory has been adopted to model the mechanical properties of several 
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BMG composites 81112,14 . It is supposed that the mechanical prop- 
erties may be changed discontinuously when the microstructural 
percolation is achieved. Although such a percolation phenomenon 
has also been suggested in the CuZr-based BMG composites in pre- 
vious studies 811 , it still lacks systematic experimental results for veri- 
fying this hypothesis solidly and the critical percolation threshold has 
rarely been achieved so far. Moreover, it still remains mysterious 
about how the microstructure forms percolation as well as how 
and why the percolation affects the tensile properties of the compo- 
sites. In this study, the tensile properties of a series of CuZr-based 
BMG composites with different microstructures have been examined 
and correlated to the microstructures systematically. The percolation 
is experimentally achieved by matching the characteristic micro - 
structural length scales, giving solid evidence for this phenomenon 
in BMG composites. The underlying mechanisms for the percolation 
and the property variations are further unraveled. It is revealed that 
the integrated tensile properties of the present BMG composites 
indicate a notable breakthrough of the general trade-off between 
strength and ductility. 

Results 

The tensile properties and fractographic morphologies have been 
systematically evaluated and analyzed for the composites with dif- 
ferent microstructures. Here we focus mainly on the structure-prop- 
erty relationships, and the readers are referred to the supplementary 
data for the detailed experimental results of representative samples. 
The dependences of the tensile yield strength o y and ultimate tensile 
strength <j uts on the crystalline volume fraction V c in the composites 
are shown in Figure 1. It is seen that o y decreases monotonously with 
the increase in V c throughout the entire crystallinity range. It has 
been revealed that the CuZr(B2) crystals precipitate from the alloys 
in a polymorphic manner that the compositions of the crystals and 
glassy matrix remain almost the same and identical to the nominal 
composition of the present alloy 815 . As a result, the properties of the 
glassy matrix can be regarded to keep constant in the BMG compo- 
sites with different V c . On the other hand, the yielding of the 
CuZr(B2) phase is supposed to correspond to the initiation of its 
martensitic transformation 1216 . Such a process depends in principle 
on the electronic structure and stacking fault energy of the crystals 
which are closely related to the compositions 16 " 18 . Thus it is reas- 
onable to postulate that the yield strength of the CuZr(B2) phase 
also keeps constant in the composites as its composition remains 
unchanged. Therefore, based on the composite mechanics, a simple 
rule-of-mixtures can be adopted to elucidate the dependence of o y on 
V c in the composites following 

Oy = O a Va + OcVc = <T«(1 - V c ) + G c Vc, (l) 

where V and a represent the volume fraction and yield strength of the 
constituents, and the subscripts a and c refer to the amorphous and 
crystalline phases, respectively. It is seen that the this approach is 
capable of describing the overall varying trend reasonably well 
(Figure 1(a)). Yet it is worth noting that there seems exist apparent 
asymmetry between the yield strengths under tension and compres- 
sion for the composites 8 . It is reported that the rule-of-mixtures fails 
to capture the compressive yield strength when V c exceeds —50% 
wherein an alternative load-bearing model 
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has been suggested 8 . However, this model is found to unreasonably 
underestimate the yield strength in tension, as indicated by the 
dashed pink line in Figure 1(a). In addition, careful examination 
further reveals that o y appears slightly higher than that predicted 
by the rule-of-mixtures in the vicinity of —50% crystallinity as desig- 
nated by the circled region. This is expected to be an indicative of the 
significant microstructural variations which will be discussed later. 
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Figure 1 | Variations of tensile strength as a function of crystalline 
volume fraction, (a) The tensile yield strength o y can be roughly captured 
by the rule-of-mixtures. (b) The ultimate tensile strength o UT s remains 
almost constant for the composites with the crystalline volume fraction V c 
lower than —60% and then decreases steeply and linearly with the further 
increase in V c . 

In comparison, g US t remains almost constant around —1700 MPa 
for the composites with V c lower than —60%, which is analogous 
to the reported case for the compressive fracture strength 10 , and 
then decreases steeply and linearly with the further increase in V c . 
The work-hardening capability can be preliminarily manifested by 
the stress window for the work-hardening process, i.e. the span 
between a UTS and a y , Act = o UTS — a y . It is seen that Act increases 
almost linearly with the increase in V c for the glass -dominated 
composites with V c not exceeding —50%. This indicates an 
improved resistance of the composites to the macroscopic unstable 
flow localization. 

Figure 2 shows the variation of the tensile ductility e T as a function 
of V c in the composites. It turns out that e T exhibits a unimodal 
hump-shaped relation with V c . It reaches its maximum value exceed- 
ing — 10% in the crystallinity range of — 40%-70% and decreases 
steeply on both sides of the plateau. This can be described by adopt- 
ing the percolation theory which is usually employed in understand- 
ing the discontinuous changes of physical properties in various 
random systems 19 . Recent studies have revealed that there exists a 
topological transition at a statistically critical microstructural con- 
dition in BMG composites where the microstructure is the most 
effective in hindering the propagation of shear bands 811,14 . This 
transition can be termed as percolation and the transition point as 
the percolation threshold. Here by treating V c as the key microstruc- 
tural parameter, e T can be modeled quantitatively by referring to the 
expressions of the percolation theory: 
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Figure 2 | Dependence of tensile ductility on crystalline volume fraction. 

The tensile ductility e T can be quantatively modeled using the percolation 
theory as denoted by the dashed curves. Representative microstructures of 
the composites with different crystalline volume fractions V c are illustrated 
in the insets. 



e T cc(vP-V c ) 



(3) 



where V p c is the percolation threshold or the critical crystalline 
volume fraction to form microstructural percolation, and /? is a 
power exponent 19 . The calculated results are proved to give a good 
description of the experimental data when the parameters V p c and /? 
are quantified as 55% and 2 respectively in the composites. Moreover, 
it is noted that the whole crystallinity range has been covered sys- 
tematically and the percolation threshold has been successfully 
achieved in the present study. These provide the solid evidence for 
the percolation phenomenon in the BMG composites. 

Discussion 

The variation in g US t can be correlated to the microstructural trans- 
ition of the composites. The composites with low crystallinities can 
be regarded as the crystalline particle-ductilized BMG composites 
with CuZr(B2) crystals dispersed within the continuous glassy 
matrix. a UST is thus dominated by the shear fracture of the matrix 
which is essentially caused by the catastrophic propagation of shear 
bands 20 . When V c exceeds —50%, in comparison, the crystals begin 
to interconnect severely to form a crystalline skeleton in the compo- 
sites and accordingly give rise to abundant grain boundaries. In this 
case, the overlapped crystals can be seen as the matrix and the glassy 
phase serves as the reinforcement. Then the global failure of the 
composites is controlled by the intergranular and transgranular frac- 
ture of the polycrystals 1011 . As a result, o UST depends not only on V c 
but also on the grain size which is also crystallinity- related 1 . 
Meanwhile, the increase in the density of grain boundaries is sup- 
posed to weaken the composites significantly. These account for the 
steep decrease in a UST with the further increase in V c . 

To uncover the origins of the percolation, the microstructural 
length scales of the crystalline particle size d and interparticle spacing 
k have been measured for the composites with crystallinities lower 
than 50%, i.e. the glass -dominated composites. It is difficult to quant- 
ify their values for those with higher V c due to the interconnection of 
individual particles. These length scales can be correlated to a char- 
acteristic material parameter of the glassy matrix-the plastic zone 
size ahead of a crack tip R p , which can be expressed in terms of the 
yield strength and fracture toughness as 



R» = K; 



ic 



(4) 



where iC IC is the plane-strain fracture toughness 4,6 . It is supposed 
that the shear bands tend to be stabilized against catastrophic 
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Figure 3 | Variations of microstructural length scales as a function of 
crystalline volume fraction. The crystalline particle size d increases 
gradually while the interparticle spacing k decreases monotonously with 
the increase in the crystalline volume fraction V c . Their ratio dlk follows a 
positive linear relation with V c , as shown in the inset. 

propagation when their mean free path for propagation is smaller 
than this critical value. For the present system, R p can be estimated to 
be -240 um by adopting K Ic ^70MPa^/m 2h22 . 

As shown in Figure 3, k decreases monotonously with the increase 
in V c and falls below R p when V c exceeds —10%, while d increases 
gradually but locates beneath R p all the way. In this manner, d and k 
approach each other as V c increases and their ratio dlk follows a 
positive linear relation with V a as shown in the inset. In respect to 
their implications, d refers to the geometrical barrier for a shear band 
to penetrate through or detour around the particles, and k represents 
the mean free path for the shear band acceleration or developing 
towards to a crack, respectively. In the vicinity of —50% crystallinity, 
the two length scales are nearly identical to each other (i.e. dlk ~ 1) 
with a value of —50 um which is smaller than the critical R p . It is 
noted that the above crystallinity conforms well to the critical thresh- 
old for the microstructural percolation. In this scenario, the majority 
of the crystals are well dispersed within the continuous glassy phase 
in the composites, resulting in an ideally inter-dispersed structure 
with the two phases partitioned by each other, as shown in the insets 
in Figure 2. On the one hand, both the two microstructural length 
scales are smaller than but on the same scale with the critical R p , 
giving rise to the predominance of stable shear banding rather than 
crack opening in the glassy matrix. On the other hand, the dimen- 
sions of the two phases match well with those of their ligaments. Thus 
the propagation of shear bands can be hindered effectively after a free 
path by the crystalline obstacles in the same scale. As a result, the 
microstructure plays the most effective role in improving the tensile 
properties of the composites. 

The effects of the ratio dlk on the mechanical properties can be 
further interpreted by simulating the stress distributions in the com- 
posites with different microstructures using the finite element mod- 
eling (FEM). Figure 4 illustrates the distributions of the equivalent 
von Mises stress in the composites with varying dlk ratios at their 
corresponding global yield strains. It turns out that the crystalline 
phase tends to deform plastically prior to the amorphous matrix 
under an applied loading due to the large discrepancy between their 
yield strengths 8 . This causes distinct misfit in the local plastic strains 
between the two phases and thus creates significant stress concen- 
trations adjacent to their interfaces 23 " 25 . It is seen that the maximum 
stress is generated near the interfaces in the direction perpendicular 
to the loading axis and extends towards the neighboring particles, 
which is analogous to the reported case in the BMGs with designed 
pores 26 . The direct results of such stress distributions can be mainly 
two-fold, i.e. to promote the initiation of shear bands at the interfaces 
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Figure 4 | Contour maps for the distributions of the equivalent von Mises stress in the composites with different microstructures. The corresponding 
global yield strains were applied for the composites with varying ratios between the crystalline particle size and interparticle spacing dlX in the FEM 
simulation. The tensile loading direction was denoted by the arrows. 



and to constrain them into the ligament regions between the neigh- 
boring particles. Meanwhile, the propagation of shear bands can be 
effectively hindered by the significant stress gradient, especially the 
steep stress drop, along the propagation path among the particles. 
The above effects on the shear bands render the composites with 
appreciable global tensile ductility in comparison with the mono- 
lithic BMGs. The cases for the distributions of the primary stress 
along the tensile axis have been further examined and are revealed 
to be almost the same or even more evident (see the supplementary 
data). 

The simultaneous increase in dlk and V c tends to enlarge the local 
strain misfit between the crystalline phase and amorphous matrix, 
and accordingly results in intensified stress concentrations near their 
interfaces. This facilitates the initiation of shear bands to a higher 
degree and makes it possible at even lower global stresses. At the 
same time, the propagation of shear bands can be hindered more 
effectively owing to the lowered global stress level and the increased 
stress gradient along their paths. As a consequence, the tensile duct- 
ility of the composites can be improved markedly as manifested by 
the left side of the plateau in Figure 2. Moreover, it is seen that such 
microstructural variation tends to decrease the angle between the 
preferred shear band plane and the global loading axis. Besides, the 
concordance between the stress concentrations on the upper and 
lower sides of each particle is gradually degraded along the inclined 
shear direction. These are also expected to retard the unstable pro- 
pagation of shear bands and prevent their premature coalescence by 
tuning their propagation path to deviate from the favored one 27 , and 
thus contribute to better ductility. With the further increase in dlk to 
exceed ~1, however, the stress concentrations around neighboring 
particles begin to overlap severely in the ligament regions. This 



results in a high stress level and decreases the stress gradient in the 
amorphous matrix and hence deteriorates the hindering effect on the 
propagation of shear bands. Such a case can also be manifested by the 
stress distribution in the composites near the global fracture strain 
which is shown in the supplementary data. Also, the ligament regions 
become too thin that can be easily penetrated by the shear bands. As a 
result, the global ductility is supposed to be deteriorated seriously in 
accordance with the right part of the plateau in Figure 2. Further- 
more, in a practical view, it is difficult to avoid the impingement of 
the crystalline particles experimentally in the rapid solidification 
process when V c is high. Thus weak grain boundaries may be intro- 
duced into the composites, which are also detrimental to the global 
ductility 11 . 

According to the above results and discussion, the CuZr-based 
BMG composites demonstrate tunable combinations of tensile prop- 
erties strongly depending on their microstructures. As the strength 
and ductility are usually major focuses in tension, a map integrating 
the two properties can be constructed for the composites. As shown 
in Figure 5(a), the composites display variable ductility e T while still 
maintaining very high strengths of o y and g US t owing to the syn- 
ergistic effects of the constituents. Their strength-ductility relation- 
ships deviate evidently from the conventional inversely proportional 
relations in engineering alloys as qualitatively denoted by the dashed 
curve 1-3 . Thus a notable exception is achieved that the mechanical 
properties can be intentionally tailored to defeat the general conflicts 
between strength and ductility and thus break through the trade-off 
between them. The pronounced combinations of high strength and 
good ductility can be readily obtained through the microstructural 
percolation by matching the length scales of crystalline particle size 
and interparticle spacing in the BMG composites. This is expected to 
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Figure 5 | Integrated tensile properties of strength and ductility for BMG 
composites, (a) The strength- ductility relationships of the present CuZr- 
based BMG composites deviate evidently from the conventional inversely 
proportional relation in engineering alloys as qualitatively denoted by the 
dashed curve, (b) The present composites are superior in tensile properties 
among various BMG composites in different systems ductilized by 
different crystalline phases 6 ' 8 " 1 1 ' 22 ' 24 ' 28 " 33 . 

provide novel strategies for defeating the general strength- ductility 
conflicts and achieving both high strength and good ductility in 
materials. Furthermore, to give an intuitive evaluation of the mech- 
anical properties, the tensile properties of other BMG composites are 
further collected from literatures and illustrated in Figure 5(b) 
regardless of their variant gauge dimensions 6 ' 8 " 11 ' 22,24 ' 28 " 33 . It is seen 
explicitly that the present composites extend a wide range of superior 
a y -e T combinations among various BMG composites. In particular, 
the combinations of variable e T with high o y , as designated by the 
dashed circle, serve as the potential candidates for tailoring the mech- 
anical properties of the composites. This usually requires the fine- 
tuning of microstructures through adopting different compositions 
or processing routes such as melt adjustment and laser surface melt- 
ing 611 ' 34 ' 35 . It is anticipated that the present composites demonstrate 
great potentials for structural applications owing to their superior 
and tunable mechanical properties and may find their practical use in 
the near future. 

In summary, the tensile strength and ductility of CuZr-based 
BMG composites show strong dependences on the microstructures 
and their microstructural dependences were described using differ- 
ent approaches. The composites exhibit notable integrated tensile 
properties through the microstructural percolation induced by 
matching the characteristic microstructural length scales of crystal- 
line particle size and interparticle spacing. The percolation theory 
has been adopted in several BMG composites, including the CuZr- 
based ones, to model their mechanical properties 8111214 , yet the 



underlying physics still remain mysterious. Here the whole crystal- 
linity range was covered systematically, providing the direct and solid 
evidences for the percolation phenomenon in BMG composites. 
Furthermore, the mechanisms of the microstructural percolation 
and the variations in the mechanical properties were unraveled in 
terms of the initiation and propagation of shear bands. The combi- 
nations of tensile strength and ductility of the present BMG compo- 
sites may help defeat the general conflicts between strength and 
ductility and break through the trade-off between them. This study 
might give deep insights on the critical issue of the structure-property 
correlations in the BMG composites and greatly aid in their micro - 
structural design and tailoring to achieve intended mechanical prop- 
erties. It may further provide new strategies for defeating the general 
strength-ductility conflicts in materials and for developing strong 
and ductile materials in the future. 

Methods 

Alloy ingots in the composition of C\x A7 £Lv A7 5 A1 5 (at.%) were obtained by arc- 
melting the mixtures of pure constituent elements under Ti-gettered high-purity 
argon atmosphere. Cylindrical rods with a diameter of 3 mm and a length of 
—45 mm were prepared by copper mold casting. To obtain various microstructures, 
the casting temperatures were varied manually by applying different melting currents 
for randomly different melting time 36 . Twenty dog-bone shaped specimens for tensile 
tests with a gauge diameter of 1.5 mm and a gauge length of 6 mm were machined 
and polished from different sections along the length of the as-cast rods. Tensile tests 
were conducted at a strain rate of 3 X 10" 4 s" 1 at room temperature in air. Here the 
tensile ductility e T was defined as the true tensile plastic strain between the yield strain 
and the total strain at the ultimate tensile strength o UTS . Morphologies of the frac- 
tured samples were analyzed using a CamScan 3400 scanning electron microscope 
(SEM). The microstructures of the specimens were characterized by examining the 
transverse sections adjacent to the same ends of the slender regions by an electron 
probe microanalyzer using the backscattered electron (BSE) mode. The volume 
fractions of the crystalline phase in the composites were estimated from its area 
proportions in the BSE images 8 " 13 . Based on the quantitative microscopy 37 , the 
interparticle spacing {X) was measured by superposing random lines on the BSE 
micrographs of the composites and the values were determined from the number of 
particles intercepting the testing line, N, and the total length between the two out- 
ermost particles (L) as k = L/(N — l) 38 . Forty random lines intercepting at least two 
particles were superposed on the micrographs for each sample. Finite element 
modeling (FEM) was conducted to simulate the stress and strain distributions in the 
composites with different microstructures. An AB-type stacking for the crystalline 
particles was considered to model the microstructures 26 . The amorphous matrix was 
treated as elastic-perfectly plastic with no work-hardening capability 23 , and a stable 
work-hardening behavior was assumed for the crystalline phase. The properties of the 
two constituent phases were collected from literatures 8 and also listed in the sup- 
plementary data. 
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